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Abstract 

Quantum phase interference and spin-parity effects are studied in biaxial 
molecular magnets in a magnetic field at an arbitrarily directed angle. The 
calculations of the ground-state tunnel splitting are performed on the basis 
of the instanton technique in the spin-coherent-state path-integral represen- 
tation, and complemented by exactly numerical diagonalization. Both the 
Wentzel-Kramers-Brillouin exponent and the preexponential factor are ob- 
tained for the entire region of the direction of the field. Our results show that 
the tunnel splitting oscillates with the field for the small field angle, while for 
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the large field angle the oscillation is completely suppressed. This distinct 
angular dependence, together with the dependence of the tunnel splitting on 
the field strengh, provide an independent test for spin-parity effects in biaxial 
molecular magnets. The analytical results for the molecular Fes magnet, are 
found to be in good areement with the numerical simulations, which suggests 
that even the molecular magnet with total spin S = 10 is large enough to be 
treated as a giant spin system. 
PACS number(s): 03.65.Bz, 75.45+j, 75.50.Xx 
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I. INTRODUCTION 



In recent years, owing mainly to the rapid advances both in new technologies of minia- 
turization and in highly sensitive SQUID magnetometry, there have been considerable the- 
oretical and experimental studies carried out on the nanometer-scale magnets which exhibit 
macroscopic quantum phenomena(MQP).0 A number of nanometer-scale particles in the 
superparamagnetic regime have been identified as candidates for the observation of MQP 
such as the tunneling of the magnetization (or the Neel vector) out of metastable poten- 
tial minimum through the classically impenetrable barrier to a stable one, i.e., macroscopic 
quantum tunneling (MQT), or, more strikingly, macroscopic quantum coherence (MQC), 
where the magnetization (or the Neel vector) coherently oscillates between energetically 
degenerate easy directions over many periods. Up to now, molecular magnets have been 
the most promising candidates to observe MQP because they have well-defined structures 
and magnetic properties. The system that have recently attracted much attention are the 
molecules [Mn 12 Oi2(CH 3 COO)i 6 (H 2 0)4] (in short Mn 12 Ac)land [(tacn) 6 Fe 8 2 (OH) 12 ] 8+ (in 
short Feg),i where tacn is a macrocyclic ligand triazacyclononance. The Mni 2 Ac molecule 
contains four Mn 4+ (S = 3/2) ions in a central tretrahedron surrounded by eight Mn 3+ 
(S = 2) ions. Oxygen bridges allow superexchange coupling among the Mn ions, and both 
high-field and ac susceptometry experiments indicate a S = 10 ground state, resulting from 
the four inner Mn 4+ spins being paralleled to each other and the other eight Mn 3+ also 
parallel with the two groups antiparallel to each other. The magnetization relaxation exper- 
iments, the dynamic susceptibility measurements, and the dynamic hysteresis experiments 
all indicate that the thermally assisted, magnetic-field-tuned resonant coherently quantum 
tunneling of magnetization occurs between spin states in a large number of identical Mni 2 Ac 
molecules.! Magnetic measurements have shown that Fes also has a spin ground state S = 10, 
which arises from competing antiferromagnetic interactions between the eight S = 5/2 iron 
spins with six spins being parallel and antiparallel to the other two spins.! 

One of the most striking effects in the magnetic MQP is that for some spin systems 
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with high symmetries, the tunneling behaviors of magnetization seem sensitive to the parity 
of total spin of the magnet. It has been theoretically demonstrated!!! that the ground- 
state tunneling level splitting is completely suppressed to zero for the half-integer total 
spin ferromagnets with biaxial crystal symmetry in the absence of an external magnetic 
field, resulting from the destructive interference of the Berry phase or the Wess-Zumino, 
Chern-Simons term in the Euclidean action between the symmetry-related tunneling paths 
connecting two classically degenerate minima. Such destructive interference effect for half- 
integer spins is known as the topological quenching.i But for the integer spins, the quantum 
interference between topologically different tunneling paths is constructive, and therefore 
the ground-state tunneling level splitting is nonzero. While such spin-parity effects are 
sometimes related to Kramers degeneracy, they typically go beyond this theorem in rather 
unexpected ways.BI~i One recent experimental method based on the Landau-Zener model 
was developed by Wernsdorfer and Sessolii to measure the very small tunnel splitting on the 
order of 10 _8 K in molecular Fe§ magnets. They observed a clear oscillation of the tunnel 
splitting as a function of the magnetic field applied along the hard anisotropy axis, which is 
direct evidence of the role of the topological spin phase (Berry phase) in the spin dynamics 
of these molecules. 

Motivated by the experiment on topological phase interference or spin-parity effects 
in the molecule Feg,i in this paper we investigate the resonant quantum tunneling of the 
magnetization vector in molecular magnets with biaxial crystal symmetry in the presence of 
an external magnetic field at an arbitrarily directed angle in the ZY plane. By applying the 
instanton technique in the spin-coherent-state path-integral represent at ion ,0 we calculate 
both the Wentzel-Kramers-Brillouin (WKB) exponent and the preexponential factor in the 
ground-state tunnel splitting. Our results show that for the small angle Oh of the applied 
magnetic field, the ground-state tunnel splitting oscillates with the field for both the integer 
and half-integer spins, and the oscillation behavior for integer spins is significantly different 
from that for half-integer spins. However, this oscillation is completely suppressed for the 
large angle region. The distinct angular dependence, together with the dependence of the 
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ground-state tunnel splitting on the strength of the external applied magnetic field, may 
provide an independent experimental test for spin-parity effects in molecular magnets. It is 
noted that the instanton technique is semiclassical in nature, i.e., valid for large spins and 
in the continuum limit. Whether the instanton technique can be applied in studying the 
spin dynamics in the molecular magnet with S = 10 (such as Fe§) is an open question. We 
study this problem with the help of exact diagonalization calculation. Our results show that 
the analytical calculation based on the instanton technique agrees excellently well with the 
exact diagonalization calculation, which strongly suggests that the molecular magnet with 
S — 10 can be treated as a giant spin system. 



II. MODEL AND METHOD 



The system of interest is a molecular magnet at a temperature well below its anisotropy 
gap, which has the following Euclidean action in the spin-coherent-state representation,^ 



S E (e,d>) = / dr 



iS (f)-is(f) cos9 + E(9, 



\drj "• \di J ' " ' ' ' (i) 

where S is the total spin of the molecular magnet. The polar angle 9 and the azimuthal 
angle <fi label the spin coherent state. E(6, <fi) is the total energy of the molecular magnet 
which includes the magnetocrystalline anisotropy energy and the Zeeman energy when an 
external magnetic field is applied. 

It is noted that the Euclidean action is written in the north-pole gauge, and the first two 
terms in Eq. flD) define the Wess-Zumino or Berry term which arises from the nonorthogonal- 
ity of spin coherent states. The Wess-Zumino term has a simple geometrical or topological 
interpretation. For a closed path, this term equals — iS times the area swept out on the 
unit sphere between the path and the north pole. The first term in Eq. ([I]) is a total 
imaginary-time derivative, which has no effect on the classical equation of motion for the 
magnetization vector, but yields the boundary contribution to the Euclidean action. Loss 
et al.EH and von Delft and Henley! studied the physical effect of this total derivative term, 



4 



and they found that this term is crucial for the quantum properties of the magnetic particle 
and makes the tunneling behaviors of integer and half-integer spins strikingly different. 

In the semiclassical limit, the dominant contribution to the Euclidean transition ampli- 
tude comes from finite action solutions of the classical equations of motion (instantons), 
which can be expressed as the following equations in the spherical coordinate system,0 

, n fd9\ . jr dE 

iS — sin0=^, 2 
\dr) d<j> V ' 

V#\ . - dE 

«D — sin = — , 3 

\drj 89' y 1 

where 9 and </> denote the classical path. Note that the Euclidean action Eq. (p]) describes 
the (1 © 1) —dimensional dynamics in the Hamiltonian formulation with canonical variables 
(f> and = S(l—cos9). 

According to the instanton technique in the spin-coherent-state path-integral represen- 
tation, the instanton's contribution to the tunneling rate T for MQT or the ground-state 
tunnel splitting A for MQC (not including the geometric phase factor generated by the 
topological term in the Euclidean action) is given by0 



T (or A) = p u p O^J 



1/2 

e" Sci , (4) 



where u p is the small-angle precession or oscillation frequency in the well, and S c i is the 
classical action or the WKB exponent which minimizes the Euclidean action of Eq. ([[]). 
The preexponential factor originates from the quantum fluctuations about the classical 
path, which can be evaluated by expanding the Euclidean action to second order in the 
small fluctuations. 

We describe the molecular magnet with biaxial crystal symmetry by the standard 
Hamiltonian,! 

n = hs 2 z + k 2 s 2 y , (5) 

where k\ > k 2 > are proportional to the anisotropy coefficients, and we take the easy, 
medium, and hard x, y, and z, respectively. If the magnetic field is applied in ZY 

plane, at an arbitrary angle ^ 9h ^ 90° with z, the Hamiltonian becomesiil 



H = hS z + k 2 S y - gfi B H z S z - gn B H y Sy, (6) 

where g is the lande factor, and /is is the Bohr magneton. The Zeeman energy term associ- 
ated with the applied field E = (0, E y , E z ) = (0, E sm9n, H cos9h) is given in the last two 
terms of the Hamiltonian. If the field is below some critical value H c {9h) (to be computed), 
the Hamiltonian Eq. (Q) has two degenerate minima, and therefore the magnetization can 
resonate between these two directions, providing a case of MQC. Under the spin-coherent- 
state and the imaginary time representation, the E(9, 0) term in the Euclidean action is 
given by 

E(9, 0) = kiS 2 cos 2 9 + k 2 S 2 sin 2 9 sin 2 — g/iBSE z cos 9 — gfisSEy sin 9 sin 

= Ki cos 2 9 + K 2 sin 2 9 sin 2 - 2K 1 (H Z /H a ) cos 9 - 2K 1 (E y /E a ) sin 9 sin 0, (7) 

with K\ = k\S 2 and K2 = k 2 S 2 being the transverse and longitudinal anisotropy coefficients 
respectively, and H a = 2Ki/gfi B S being the anisotropy field. Introducing A = K 2 /Kx, 
cos#o = E cos9h/ E a , and sin0 o sin0o = E sin 9h / ^H a , the E(9,(f>) reduces to 

E(9, 0) = K\ (cos# — cos^o) 2 + K 2 (sin 9 sin — sin 9q sin O ) 2 + Eq, (8) 

where Eq is a constant that makes E(9, 0) zero at the initial state. It is clearly shown in Eq. 
(§) that the energy minima of system are at 9\ = 9q, 0i = 0o and 9 2 = 9q, <p 2 = tt — 0o, and 
therefore there are two different instanton paths of opposite windings around hard anisotropy 
axis. We denote them by instanton A : = O — > = +tt/2 — > = 7r — O , and instanton 

B : = 0o ► = — 7r/2 ► = 7T — 0q. 

The critical field at which the energy barrier disappears can be determined by, 

Ecos9 H 

cos6» = ^ 1, (9) 

Esm9u , , 

sin <^o = — (10) 



Taking into account Eqs. (P) and ([TDj) , we obtain 



For the special cases Oh = 0, ir/2, we have H c = H a , XH a , respectively. The dependence of 
H c /H a on H is plotted in Fig. 1 for A = 0.71. 

Now we investigate the tunneling behaviors of magnetization by applying the instanton 
technique in the spin-coherent-state path-integral representation. First of all, we must find 
out the classical path and (f) that satisfies the boundary condition. Along the classical 
path, E(0, 0) is conserved, so that the relation between 0{r) and 0(r) can be found purely 
by using energy conservation.!'! After some algebra, we obtain 

(cos 0o — A 1//2 y sin <fi) ± iX 1 ^ 2 (x sin <fi — sin 0$ sin 0o) 

cosO = -, (12) 

1-Asin 2 v ; 

where x, y are variables defined by 

*={ (a2 + T + f 2 and ^{ < a2 + " 2 > 1/2 - a } 1/2 , (13) 

with a = 1 — cos 2 O — A sin 2 + A sin 2 sin 2 O and b = 2A 1//2 sin cos sin O . Here we 
assume that the condition a ^ is always fulfilled for the small magnetic field H. The 
positive and negative sign in Eq. ( |T2"D are corresponding to instanton and anti-instanton 
solutions for path A or B, respectively. For the low magnetic field , one can take 0(r) to 
be entirely real (see appendix A for detail) .i Then through the expression S c i = iS — 
cos(0))(d<f)/dT)dT, the desired Wentzel-Kramers-Brillouin (WKB) exponents (or classical 
actions) are found to be 

25'cos#o 

rvr/2 



lm{S*) = S(+n - 20 o ) - ( f_™ffi 2 (f - arctan ((1 - A^tan^)) 



+ 2A 1 / 2 ^ r 

J d>„ 1 - A sm 



2 JL> 



Re(^) = 2X^S { ^^^t (-f + arctan ((1 - A)V 2 tan, 



7r ' 2 x sin 



+ / (14) 



4>» 



for instanton A, and 
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Im(Sf) = S(-7T - 20„) + 2 1 5 _ C ° S , / ° 2 (| + arctan ((1 - A) 1 / 2 tanfo)) 

1 - A sm 



Re(^) = 2X^S { ^^p (+| + arctan ((1 - A)^tan0 o )) 
r/ 2 x sin <f)d<p \ 

+ L r^w0/ (15) 

for instanton B. Therefore, the difference of the action of two tunneling paths is 

Im(A^) = Im(S*) - Im(^) = -2Svr (l - ^y/i ) = ~ 2 ^ H )> ( 16 ) 
Re(AS d ) = Re(S%) - Re(^) = 2Sn sin 6 sin O ( j^J • (17) 



Note that Im(AS c /) comes from the Berry phase iS(l — cos 6) (dip/dr), and leads to the 
oscillation of the ground-state tunnel splitting with the magnetic field. 

The preexponential factor in the tunnel splitting can be evaluated by taking the asymp- 
totic form of zero-mode (dcp/dr) ,0 which can be deduced from the classical equations of 
motion (see appendix A). 

Finally, the ground-state tunnel splitting is found to be 



£A = frA |l + e 2 ^e- ReA ^| 

= hA { (1 - e - RcAS «) 2 + 4e" ReA ^ cos 2 ^H)}^ , (18) 



with 

° ° (sin 2 O - A) 1/2 V 2*K J Un 2 £ o + Acos^osin 2 0j ' 1 j 

where u = (2KJS) A 1 / 2 sin 9 cos 0o- The dimensionless prefactor c can often be of order 1 
or so, and is therefore relevant to the exact diagonalization calculation. 
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III. RESULTS AND DISCUSSION 



Before we discuss the Eq. flIB|) , we note here that our model can be directly related to 
the model describing the molecular Fe§ magnet,! 

H = -DS 2 Z + E(S 2 - S 2 ) + gfisH'Jx + 9^ B H y S y , (20) 

with Ki — (D + E)S 2 , K 2 = (D - E)S 2 and H z = H' x , H y = -H' y . According to the 
typical parameters of Fes, D = 0.275K, E = 0.046K and g = 2, we obtain that K\ = 32. IK, 
K 2 = 22. 9K, A = 0.71 and H a = 4.77T. These parameters are used throughout the whole 
calculation. 

First, we discuss the effects of arbitrarily directed field on the spin phase oscillation. 
Our results show that the topological phase interference or spin-parity effects depend on the 
direction of the magnetic field significantly. From Eq fll8D , whether the ground-state tunnel 
splitting oscillates with the field is determined by two factors, §{H) and Re(AS , c ^). For a 
fixed field strength, i.e. H = Hs = {\JH'y + Hfj , it brings two- fold effects to increase the 
field angle Oh'- (i) the dependence of &{H) on the field strength is reduced. As 9h increases 
up to |, <&(H) = nS (l — j^^j2ff~\ gradually reduces to a constant ttS. (m) The degree of 
oscillation a which is defined as a ratio of the oscillation part 4 e - ReAS, c; ^ Q ^he non-oscillation 
part (l — e - ReAs c^ 2 in Eq. flTSP decreases from +oo to a small value. 

For 6h = (the field is along the hard anisotropy axis), the symmetry of two classical 
paths imposes Re(AS , c ;) = 0, then, the tunnel splitting A = 2A |cos ( l ) (iJ)| oscillates with 
the field and thus vanish whenever 

| = (1 _ A)1/2 (^V2)_ (21) 

where n — 0, 1,2, ... is an integer.! For Oh — f (the field is along the medium anisotropy 
axis), on the other hand, §{H) becomes a constant ttS, and therefore the tunnel splitting 
increases monotonically with the field. When the field is applied in ZY plane with an 
arbitrarily angle Oh, one may expect to observe a crossover for a certain field angle C H . 
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Choosing the magnitude of H to be the first value that makes cos$(iJ) = and taking 
a = 1, we obtain 

9 C H « arctan ^ L7 ^ A ' \ (22) 

For the molecular Fes magnet, A = 0.71, we find 8 C H ~ 25°, which agrees well with the exact 
diagonalization calculation (shown in Fig. 2a) and the experimental results obtained by 
Wernsdorfer and Sessoli.@ It is noted that the value of 9 C H is independent of the total spin 
S of the molecular magnet, and depends on the parameter A only. For highly anisotropic 
materials, the typical values of the transverse and longitudinal anisotropy coefficients are 
K\ « 10 7 erg/cm 3 and K% « 10 5 erg/cm 3 . Thus, A = 0.01, the value of 9 C H is estimated to be 
3.2° and is much smaller than that for the molecular Fe§ magnet, which means that even a 
very small misalignment of the field with the hard anisotropy axis can completely destroy 
the oscillation. Therefore, the molecular Fes magnet is a better candidate for observing 
the oscillation of the ground-state tunnel splitting with the field compared with the highly 
anisotropic materials. 

Next, we turn to the comparison between the analytical calculations and the numerical 
stimulations. It is noted that the instanton approach is semiclassical in nature, i.e., valid 
for large spins and in the continuum limit. Whether the instanton technique can be applied 
in studying the spin dynamics in molecular magnets with S = 10 (such as Fe§) is an open 
question. We have performed the numerical diagonalization of the Hamiltonian Eq. (^|) 
for the molecular Fe§ magnet in the presence of an external magnetic field at an arbitrarily 
directed angle in ZY plane. As illustrated in Fig. 2a, for Fe 8 , the analytical results based 
on the instanton technique are in good agreement with the exact diagonalization results. In 
order to show the agreement more clearly, we have presented the numerical and analytical 
results in Table 1 for the molecular Fe§ magnet in the magnetic field applied at angle 
9jj = 0°, 15°, 30° and 90°. It is clearly shown that the accuracy of the semiclassical 
calculation is very high for the low magnetic field up to H = LOT for the entire region 
of the angle 0° ^ 9h ^ 90°. As a result, we conclude that the molecular magnet with 
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S = 10 can be treated as a giant spin system. The numerical and analytical calculated 
tunnel splitting as a function of the field are shown in Fig. 2b for half-integer spin S = 9.5, 
and the good agreement between numerical and analytical results is also found. From Fig. 
2, it is obviously that the tunneling behavior of magnetization of integer spins is significantly 
different from that for half-integer spins. At the end of this section, we present the results 
with different parameter A for the field along z or y in Figs. (3a) and (3b). It is interesting 
to note that the tunnel splitting increases rapidly by lowing the parameter A, which suggests 
that highly anisotropic materials are more suitable for observing MQP in experiments. 

In conclusion, we have studied the ground-state tunnel splitting in the molecular magnets 
with biaxial crystal symmetry in the presence of an external magnetic field at an arbitrarily 
directed angle. The switching from oscillation to the monotonic growth of the ground-state 
tunnel splitting on the field angle has been shown in detail. Our results are suitable for a 
quantitative description of some aspects of the new experimental behavior on the molecular 
Fe 8 magnets.i 



The authors are indebted to Professor W. Wernsdorfer and Professor R. Sessoli for pro- 
viding their paper (Ref. 9). The financial support from NSF-China (Grant No. 19974019) 
and China's "973" program is gratefully acknowledged. 



In this appendix, the general scheme for calculation of the ground-state tunnel splitting 
is presented, the main assumptions and approximations are also outlined. 

We start with the relation between 6{r) and 0(r) obtained from the energy conversion^ 
(see Eq. (0)). Only one instanton, say the + one, need be considered explicitly (from now 
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Appendix A: Evaluate the ground-state tunnel splitting 



on, we drop all the bars and identify 9 = 6{r) , <p = <K r ) f° r convenient): 




(23) 



1 — A sin 2 (f> 
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(x — A sin 9 sin O sin 0) — i (A 1//2 cos 9 sin — y) 

sinfl = , 2 V . (24) 

1 — A sm 

In order to determine the classical paths, we need another equation of motion(see Eq. @), 
dE(6,4>) 



iS sin i 



86 

2Ki(cos6 — cos^q) sin — 2K 2 (sm6 sin — sin 8q sin O ) cos # sin 0. (25) 



After dividing the factor 2J^ 1 sin6 ) on both sides and substituting Eqs. ( Effi) and (0), we 
rewrite Eq. ( |25"D as 

. 5 , 



—\ l l 2 y sin + iA 1//2 (x sin — sin Q sin O ) + 
A sin 9 Q sin O sin < 



'IK?' 

(cos #o — A 1 / 2 ?/ sin 0) + zA 1//2 (x sin — sin #0 sin 1 



(x — A sin # sin O sin 0) — % (A 1 / 2 cos 9 sin — y) 
Re/(0)+2lm/(0), (26) 



where 



Re/(0) = — A 1 ^ 2 ysin0 + A sin #0 sin O sin 

(l — A sin 2 0) [x cos 80 — A 1 / 2 sin 6* sin O ) 
(x — A sin 8 sin O sin 0) 2 + (A 1 / 2 cos #0 sin — y) 2 



Im /(0) = A 1 / 2 (x sin — sin # sin 0o) + A sin # sin O sin 

[(cos #0 — ^^y sin 0) (A 1 / 2 cos # sin — y) + A 1 / 2 x 

(x sin — sin 9q sin O ) (x — A sin #0 sin O sin 0)] 

x — . (28) 

(x — A sin #0 sin O sin 0) + (A 1 / 2 cos 6*o sin <p — y) 

As shown in Eq. (^) , if Re /(0 ) is always zero for arbitrary 0, one can then take 0(r) to 

be entirely real. Unfortunately, it is not the case, Re/(0 ) is not always zero except H z = 

or H y = 0. Checking the right-hand-side of Eq. ( f27f ) more carefully, we find that up to the 

leading term, Eq. (^) can be rewritten as 



10 f / 1 \ 1 A 2 sin 3 . 

Re/(0) «- . cos flo am 

1-Asin 2 ^ 7 
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A 2 

,1/2 



^ -TTj cos O sin cp 

(1-A) ] 



A H z H y 



;i-a) 1/2 hi 

A H 2 sin 20 H 



(29) 



2 (1 - A) V2 #a 

It is clearly shown from Eq. (|29"D that the Re /(0) is almost two orders of magnitude smaller 
than Im/(0) for the low magnetic field(i.e., for H/H a ^ 1/5). So it is sufficient to drop 
the term Re /(0) in Eq. (^) for the first order approximation and treat 0(r) as a real 
variable. Under this assumption, the classical Euclidean action can be evaluated easily 
without resolving the analytical solutions of 0(t) and 0(t), 



S d = iS (l-cos(0(r)))(d(j)/dT)dT 

J — oo 

= iS [^(1 - cos(0(4>)))d<t> 



r+f ( (cos 6»o - X 1/2 y sin 0) + zA 1/2 (x sin - sin O sin O ) , , 
= iS / 1 - . 2 (30) 

where 0j = 0(r — > — oo), 0y = 0(r — > +oo). The results for instantons A and B are given 
in Eqs. <$L§ and flTJ). 

Next let us study the preexponential factor at two different magnetic field directions of 
H = O and tt/2 ^ H > 0. 

I.0 H = O 

The preexponential factor can be evaluated from zero-mode (^/dr)fl It is obvious from 
symmetry that the preexponential factors along two instanton path A and B are equal. 

For H = 0, we have O = and 0, = 0, 0/ = tt. Then, from Eq. (|13|) we get 
x = (l — cos 2 O — A sin 2 0) and ?/ = 0. Therefore, Eq. (|26|) reduces to 



= A l/2 ^ _ cog 2 0q _ X gin 2 ^ V2 gin _ (31) 

This equation can be integrated easily. Defining to = 2 -| J -A 1//2 sin^o, we obtain 
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coso= (l — cos 2 9q — A) — _ , ° , 9/ ^ 1/2 - (32) 



1/2 tanh(u; r) 

(l — cos 2 #o — A tanh 2 (wot)) 

It is easily verified that 0^0, 7T, asr^ioo. 

Following the standard procedure of the Ref. 10, we write the final result as 

hA = cwf 85ing ° 1/9 ( JL.) 1/3 e -B4 \ C0S ^ H )\ , (33) 

(i_ cos 2^ _ A) i/2 \2ttKJ 1 ^ 



where 



(A 1 / 2 \ / 1 -L- A 1 / 2 cos 6>o 

1± ff) — ^.lu (i :r s r° i > (34) 
1 _ A 7 / Oh _ xW 2 I 1 - A 1 / 2 cos go I ' 
1 sin 0o / Z ^ A ^ V 1 (l-A^sine,, 



^H)=nS[l--^-\. (30) 



and 



(1 - A) 1 / 2 

Eq. (|3~3"D can be also deduced from Eqs. (0), flT6|) , ( |17D and (pT8|) for the special case O — 0. 
It is interesting to note that Eqs. (0) and ( |3~5| ) agree exactly with the previous result (Eqs. 
(3.10) and (3.11) in Ref. 7) found by Garg for the molecular Fe§ magnet. 

II. vr/2 > 9 H > 

For n / 2 ^ > 0, it is difficult to integrate the equation of motion Eq. (|26|). So we have 
to take the asymptotic form of zero-mode (dip/dr). Notice that — > tt — 0o as r — > +oo, we 
can expand Im/(0) according to the small parameter a = <p — (tt — O ). Up to the leading 
term, Eq. fl26|) has the form 

= -w (0 - (vr - O )) , (36) 

where wq = 2 -| 1 A 1//2 sin^o cos0o- We then integrate the Eq. ( p6|) 

= vr - 0o - c e'^ T , (37) 

where Co is an integration constant. Therefore, the asymptotic form of zero-mode (dtp/dr) 
can be written as 
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t = c uj e 



-U> Q T 



(38) 



Straightforward, following Ref. 10, the tunnel splitting hA A and hA B corresponding to path 
A and B have the form 

/ S 2 \ 1/2 ( sin 2 9, 



CaU ° V 27r ^J Vsm 2 ^o + A cos 2 sin 2 0o 



hA = cbuJq 



3/2 



s 2 



1/2 



sin 2 9 



2ttKi J y sin 2 9o + X cos 2 #o sin 2 0o 



1/2 



Re(S c 4) e -iIm(S c 4) 



-Re(5f) -ilm(S^) 



(39) 



(40) 



where ca and are the numerical factors in order of 0(1) that come from the integration 
constant mentioned above (cq in Eq. fl38f) ). Because the existence of H y breaks the symmetry 
between the classical path A and B, ca may not be equal to cb- But for the small magnetic 
field, we assume theirs equivalence. Comparing Eqs. (|39| ) and (pUf) with Eq. (p3|), we find 

4 sin 9 



c A = c B = c- 



(41) 



(l-cos 2 e -\) 1/2 ' 

Now we turn to derive the total tunnel splitting hA using a recently proposed effective 
Hamiltonian approach. For the present case, the effective Hamiltonian is found to be 



H 



eff 



(hA A + hA B ) 

(hA A + hA B )* 



(42) 



Diagonalizing the effective Hamiltonian, we obtain the desired result as shown in Eqs. (18 
and (jlj). 
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e H 


= 0° 


0h = 


15° 


0H = 


30° 


H = 


90° 






A d 


A, 


A d 


A, 


A d 


A, 


A d 


A, 


H s 


= O.OT 


0.00683 


0.00683 


0.00683 


0.00683 


0.00683 


0.00683 


0.00683 


0.00683 


H s 


= 0.2T 


0.00556 


0.00558 


0.00778 


0.00780 


0.01445 


0.01448 


0.05651 


0.05658 


H s 


= 0.4T 


0.00165 


0.00167 


0.01714 


0.01733 


0.06521 


0.06591 


0.75350 


0.75752 


H s 


= 0.6T 


0.00502 


0.00515 


0.04658 


0.04782 


0.29419 


0.30232 


8.20514 


8.30553 


H s 


= 0.8T 


0.01417 


0.01485 


0.13101 


0.13783 


1.34388 


1.42229 


73.1429 


74.7913 


H s 


= LOT 


0.02362 


0.02548 


0.39753 


0.43399 


6.19230 


6.88449 


537.366 


557.203 


H s 


= 1.2T 


0.02505 


0.02808 


1.31263 


1.51459 


28.5282 


34.3140 


3277.23 


3462.35 


H s 


= 1.4T 


0.00799 


0.00942 


4.70886 


5.90363 


129.907 


176.358 


16704.7 


18092.3 



Table 1. The ground-state tunnel splitting A d (calculated by diagonalization) and A; 
(calculated by instanton approach) with 9 H = 0°, 15°, 30° and 90° are listed as a function 
of magnetic field Hs- The units for the tunnel splitting and magnetic field are 10~ 7 Kelvin 
and Tesla, respectively. Here, S = 10, K x = 32. IK, K 2 = 22.9K and A = 0.71 are used for 

the molecular Fe 8 magnet. 
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Figures Captions 



Fig. 1. The critical magnetic field H c is plotted as a function of anlge 6 H at A = 0.71. 

Fig. 2. The ground-state tunnel splitting A are plotted as a function of the magnetic 
field H s at H = 0°, 3°, 10°, 30°, 60°, and 90° for (a) S = 10 and (b) S = 9.5. The 
other parameters are K 1 = 32. IK, K 2 = 22. 9K and A = 0.71, the typical values of the 
molecular Fes magnet. The results of the instanton approach and the exact diagonalization 
are represented by the solid lines and square symbols, respectively. 

Fig. 3. The ground-state tunnel splitting A are plotted as a function of the magnetic 
field H s for A = 0.3, 0.5 and 0.7 with S = 10, K x = 32.1K at (a) 6 H = 0° and (b) 9 H = 90°. 
The results of the instanton approach and the exact diagonalization are represented by the 
solid lines and square symbols, respectively. 



18 



